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fisclude <stdio.h> 

void 

nainO 

{ 

float 

dC4] [3] , a[3] [3] , K3] , cC3] , ddC3] .r[3][3] ,kkC3] [3] , 33, sss; 
int i.j.k; 

A input X & y of VSSi */ 

printfClSPDT RGB and Hhitrfn"); 

printfClx Ry Gx fiy to By Ix 1fy%i".); 

scanfCSf Sf Sf Xf Xf Xf Xf Xr, &dt0l[0],4d[0][l] 

. ad[l]t03,M[l]Cl] 
, W[2][0],4d[2][l] 
. W[3]C0],U[3]Cl]); 

/* 



d[O][0] 




0.67; 




d[0][l] 




0.33; 




d[l][0] 




0.21, 




dtllCl] 




o.n, 




d[2][0] 




t).14, 




d[2][l] 




0.08, 




dC3]C0] 




0.31, 




dC3][lI 




0.316; 



«/ 



FIG. 20 



/* calculate z fraa z & 7 «/ 
for(i = 0; i < 4: i++){ 

if((d[i]CO] + d[i][l]) > l.OH 
d[i](2] = 0.0; 

} 

d[i][2] = 1.0 - d[i][0] - d[i][ll: 



priii«(".MAIRIX&'); 
ford = 0; i < 3; 

printfC^f ); 

for( J = 0: j < 3; 

printf(-X5.3f?f .dCi][j]): 

} 

printfrYn*); 

} 
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/« ealacnlate utriz «/ 
{ 

int 11, 12, jl, i2; 
ford = 0: i < 3; 

11 = i + 1; 

12 = i + 2; 

If (11 > 2) il = 0; 
if (12 > 2) 12 = 12 - 3; 
for(j = 0; j < 3; i++K 
il = j + 1: 
j2 = j + 2: 
if (jl > 2) jl = 0; 
if (j2 > 2) j2 = j2 - 3; 

aCnti] = d[il][jl]»dCl2][j2] - dCil][j2]«dCi2][jl] 

} 

} 

} 

/* calenlate of BUNBO «/ 
ford = 0;i < 3; l++){ 
b(l] = 0; 

for{j = 0; j <3; j++){ 
bCi] = a[i]Cjl . d[3][j] + bCi]; 

} 

} 
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/« MAHIX »/ 
forU = 0; i < 3; 
for(j = 0: j < 3; j++){ 
a[i][j] = a[i][j] / b[i]; 
rCi][j] = a[i][j]; 
if(i = iH 

UtCi][j] = 1.0; 
} else { 

kkCi]tj] = 0.0; 

} 

} 

} 

/« IH7ESSZ HAIEIZ «/ 
for(i = 0; i < 3; i++H 
for(j = 0; j < 3; j-H.){ 

ddCil = atj][i]: 

a[j][i] = 0.0; 

} 

a[il[i] = 1.0: 
for(j = 0; i < 3; j-H-){ 
c[j] = a[i][j] / dd[i]; 

} 

forU = 0; j < 3; i++){ 
for(k = 0; k < 3; !c++){ 
aCj][kl = a[j][k] - c[k]sdd[j]; 

• } 
} 

for(j = 0; j < 3; i++K 
a[i][j] = c[i]; 

) 

} 
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/« SEinKA «/ 

ss = a[l][0] + aCnCl] + aCl][2]: 
sss = rCl][0] + r[lKl] * rCl][2]; 
for(i = 0; i < 3; i++){ 

for(j = 0; j < 3; j++}{ 

a[i]tj] = a[i][j] / ss; 

rCi][j] = rCiKj] / sss; 

} 

} 



FIG. 2A 



/* resnlt «/ 

printf ("original dat^'); 
for(i = 0; i < 4; i«-){ 

printf("?r); 

for( a = 0: j < 3; j++)-{ 

printfCXr.Sf ",d[i][j]); 

) 

printfCfti'); 

} 

printfCHAHUHii'); 
foKi = 0; i < 3; i++){ 

printf(-?tr); 

for( j = 0; j < 3; j++){ 

priatf(-X7.5f ".rCilti]): 

} 

printfr?!!"); 

> 

printf r IHVIRCE HATaHfti'' ) : 
. foKi = 0; i < 3; i+4){ 

printf(-¥t-); 

foK j = 0; j < 3; j++){ 

printfrrr.Sf -,a[il[jl): 

} 

printf (-in*); 

} 

for(i = 0; i < 3; i++){ 
for(j = 0: j-< 3; j++H 
kk[i][j] = *CiltO]*rCO][i] + a[il[l]^[l][i] + iCi][2]*r(2]Cj3; 

} 

} 

priatf(-IAmAIf!n-); 
foKi = 0; i < 3; i-H-H 

printf (-¥t"); 

for( j = 0; j < 3; j++){ 

printfrrr-Sf ",kk[i][j]): 

} 

printfrVn"),- 

} 



